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УДК 33:519.2 
А.Ю.ЛУКОНИН, Ю.А.ЛУКОНИН 


СТОХАСТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ПРОДАЖ 
С ПОМОЩЬЮ СЛОЖНЫХ ПУАССОНОВСКИХ ПОТОКОВ 


Случайный процесс продаж рассматривается как случайный пуассоновский поток 
коммерческих сделок, размеры которых также являются случайной величиной. По- 
казано, что если известно распределение вероятностей размеров сделок и плот- 
ность потока покупателей, то можно вычислить распределение вероятностей раз- 
меров продаж. Сформулированы условия, при которых это распределение продаж 
не будет нормальным. Приведены примеры негауссовских распределений в реаль- 
ных процессах продаж. 

Ключевые слова: случайные процессы, случайные пуассоновские потоки, распре- 
Деление вероятностей размеров продаж, микроэкономика фирмы. 


Введение. Процессы продаж являются разновидностью потоков коммер- 
ческих сделок, в которых в основу стохастических моделей положен 
Пуассоновский характер потока сделок, характеризующийся как случайный 
нестационарный поток независимых событий без последействия. 


Обозначим 1 (1) - плотность потока покупателей. Тогда количество 
продаж 7 за интервал времени длиной {, начинающийся в точке &, под- 
чиняется закону Пуассона 


т 


Р.(т,1)= — (т= 0,1,...), (1) 
т. 


где 4 - математическое ожидание числа продаж на участке [1,4+1], 
+1 
а= 1(04. (2) 
Ио 
Распределение вероятностей р, потребности покупателя в А еди- 


ницах товара будем называть в дальнейшем распределением потребностей 
покупателя в товаре. Тогда случайная величина /Л, равная количеству 


проданного товара за интервал [1,2 +7], может быть выражена следую- 
щим образом: 


М№ = К, ' (3) 


где Т - случайное число покупателей на отрезке [1,1 + 7 ] , имеющее рас- 


пределение (1), а А, - количество товара, купленное /-м покупа- 
телем: случайная величина, имеющая распределение вероятностей 
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р,., К=12,... (Ррь - вероятность того события, что один покупатель 
купит Кединиц товара). 
Можно определить математическое ожидание и дисперсию случайной ве- 
личины М [1]: 
= * 
М[М]= М[]*а, (4) 


ЫМм]= ((М[К])” + БА] *а , (5) 
где Мр [*] - математическое ожидание и дисперсия соответствен- 
но; величина @ задается выражением (2); 


МЕ К*р; Ма МУ * 
й=1 К=1 

Соотношения (4) и (5) являются достаточными для описания стати- 
стических характеристик сложного пуассоновского потока Л’ [1] только в 
том случае, когда величина /№ распределена по нормальному закону. 
Именно на предположении о том, что величина № распределена нормаль- 
но, строится современная логистика управления запасами. Однако в рабо- 
те [2] это фундаментальное предположение подвергается сомнению и при- 
водятся примеры, когда эмпирические оценки распределений вероятностей 
сложных товарных потоков мало похожи на гауссовские распределения. 

В работе [1] отмечается, что сложные пуассоновские потоки имеют 
сложные распределения, не имеющие аналитического выражения. Тем не 
менее, такие распределения нужно рассчитывать для построения систем 
управления товарными и финансовыми потоками фирмы. Поэтому рассмот- 
рим алгоритм вычисления функции распределения вероятностей для слу- 
чайной величины №. 

Распределение вероятностей для сложных пуассоновских пото- 


ков. Пусть распределение потребностей покупателя Р, определено на 


всей числовой оси - К С:РИриК ий, ПрШЕ > 0, а 


поток покупателей имеет характеристики (1) - (2). Тогда вероятность 
того, что 7 покупателей купят Л’ единиц товара, равна 


О, (М) = Рь Рь--- Ре Ру-ь-ь-- (6) 


в =- №=- Аи-1=- 
а вероятность того, что на интервале [%,+ ТГ] будет продано № единиц 
товара, равна 
РМ;к,Т)= Р‚ (Тк )О„(М), (7) 
т=1 


где Р.(Т,&,) - вероятность появления Т покупателей товара на интерва- 


ле [4,к + 7], определяемая соотношениями (1) - (2). 
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В правой части выражения (6) сделаем замену переменных: 
а 
о 
а. 
И аа (8) 
А Е аа 


Вы = А = И = 1,2 в а я 
Тогда выражение (6) можно переписать в виде 


О, (М) = Рм-ь Ри-ь--. Ра, ,-1 


т-1 


Р‚_. (9) 


1=- Ь=- и-1=- 
Учитывая, что р, = О при К 0, нетрудно видеть, что в (9) будут отлич- 
ны от 0 только те слагаемые, для которых выполняется система неравенств 


20: 
е = т > 0; 
НИ 22-0: 
т-3 2 (10) 
И: 0: 
М-1> 0. 
Нетрудно видеть, что необходимым условием, при котором удовлетворяет- 
ся система неравенств (10), является положительность 
11= 1,2,...тин № . Тогда (9) можно записать в виде 
О (М)= Ре” Ре ПИ М> 0. (11) 
1=1 Ь=1 1.1=1 
В выражении (11) перейдем к матричной записи. Введем столбцы 
0}, =О, (1), 1= 1,2,.... (12) 
и 
Рева (13) 
и матрицу 
вы р. уе (14) 
Тогда (11) запишется в виде 
т-1 
(0,}=[Р]” {в}. (15) 


Из области определения р, следует, что [Р р, = 0,Г0ег 7, тогда в выра- 


жении (15) матрица [Р] имеет форму нижней треугольной матрицы с ну- 


левой диагональю. Рассмотрим случай, когда множество значений [, В КО- 
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тором р,> 0, конечно и ограничено сверху / Г (ядро Р). Тогда нижний 


треугольник матрицы [Р] вырождается в ленту шириной Г, и матрица 
имеет следующий вид: 
0... 


ро... 
р» 210... 


М 5 


Е р 0... 
О р, ре, 
0:0 РЕ р 0 


Матрица Р обладает следующим свойством: результатом умноже- 
ния слева матрицы Р’ на вектор с границами ядра // / Г, будет также 


вектор с конечным ядром, имеющим границы 2+1 / Г, + Г. Таким об- 


разом, в результате воздействия матрицы Р на вектор с конечным ядром 
ядро растягивается и сдвигается на 1. Если распределение вероятностей 
сложного —пуассоновского потока представить в форме столбца 


[Ра,Т)}, = Рец, Т),1= 1,2,..., тогда из (1),(7),(12) и (15) следует, что 


Рот - ео Ц. а 


Численная оценка распределения вероятностей. Для того чтобы оце- 
нить количество слагаемых М, при котором конечная сумма в правой ча- 


сти (17) с заданной точностью # аппроксимирует [Ра’Т) ‚ оценим нор- 
му разности векторов 


Вт" р 

















т=1 т! 
(„Ту ее 
_ р-аТ) а(, р т-1 
` т=М+1 т! | [и 
Очевидно, что 
КР = [21 (19) 
1=1 
и 
т 
[2] = пах [Р],= Р=1. (20) 
Л 
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Отсюда 
(а@,Т))” 


т=М+1 т! 


(21) 








ЧИ] 


т= М+1 т! 





Правая часть в (21) есть не что иное, как остаточный член А/(х) разло- 


жения е” в ряд Маклорена. Представим А„(х) в форме Лагранжа [3]: 
М+1 


Ре 8 
Ки(х)= Пил а 0‹9 1. (22) 


Тогда для правой части (21) получаем оценку 
(а(к,Т))" (а@,Т))“ рат) 
т= М+1 т! (МП! 
Обозначим и= М +1. Для оценки п! используем формулу Стирлинга [3] 


(23) 


п 


9 
п!= 21 п ф еёде , Ох < 1, (24) 


е 
откуда получаем оценку 
(а%,Т)) (а&,Т)у 
(п)! п’. (25) 
211 — 
е 
Используя (18), (21), (23) и (25), получаем, что если найти такое ми- 
нимальное значение И , при котором выполняется неравенство 
а(,Т)е ь 
п 
тогда будет также выполнено неравенство 


[Рарт}- е®” Е 


е\/2и, (26) 


# : 127 











т=1 
Обсуждение результатов. Используя метод половинного деления, для 
заданного ё можно легко найти ши(и), при котором выполняется нера- 


венство (26), в зависимости от а(&,Г). Например, если обозначить 

М = [а(%,Т)*е]+ АМ (а] —целая часть а), то при & = 0,001, можно по- 

лучить зависимость д М от а(®,Т)*е , представленную в таблице. 
Зависимость А М’ от плотности потока покупателей 


2 382 


383 2909 
2910_21537 
21538 ОИСИ 
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Таким образом, видно, что поправка Л М существенна лишь при малых 
значениях а(к,Т)*е 

При оценке распределения продаж необходимо учитывать, что со- 
гласно центральной предельной теореме при а(®,Г)> оно должно 


стремиться к нормальному. На деле если ядро распределения потребно- 
стей покупателя локализовано в пределах 0 20, то уже при значениях 


а(к,Т): 7 19 вид распределения продаж практически не отличается от 


нормального. 

Если же распределение потребностей покупателя в товаре имеет 
«длинный хвост», то, например, распределение величин ежедневных про- 
даж может существенно отличаться от гауссовского. На рис.1 показано та- 
кое распределение некоторого товара сантехнической группы, полученное 
по выборке из 2346 покупок, предоставленной торговой компанией ЗАО 
«Билд». 











оо ры ®) юятооон 
и 











Рис.1. Распределение потребностей покупателя 


Ядро этого распределения локализовано в интервале 0< Г, 420. На ри- 
сунке хорошо видны пики распределения, соответствующие покупкам, раз- 
меры которых кратны 7. Это объясняется тем, что этот товар в основном 
покупали упаковками по 7 штук. Ряды покупателей этого товара и продаж 
за 2004 год показаны на рис. 2 и 3. На них также представлены тренды 
этих рядов, полученные методом «скользящего среднего» с шириной окна, 
равной 30 дней. На рисунках хорошо виден нестационарный характер этих 
рядов. Наличие отрицательных значений в ряде продаж объясняется тем, 
что в данном представлении учтены возвраты товара покупателями, кото- 
рые учитывались с отрицательным знаком. 





о} 








я 





Рис.2. Временной ряд ежедневных продаж 
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Рис.3. Временной ряд покупателей 


На рис.4 и 5 представлены гистограммы распределения вероятно- 
стей ежедневных продаж 07.05.04 и 22.09.04 соответственно, вычисленные 
с погрешностью 0,001. Оценки плотностей потока покупателей в эти дни 
соответственно 5,8 и 15,07 покупателей в день. 
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Рис.4. Распределение продаж от 07.05.04 
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Рис.5. Распределение продаж от 22.09.04 


На рис.4 хорошо видно влияние упаковок на распределение про- 
даж, на рис.5 это влияние уже нивелировано. Расчеты проводились в си- 
стеме «1С предприятие» на РС РепЧит-Т\ 1,2 Ггц, время расчета 2 мин и 
3 мин 34 сек соответственно. 
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Выводы. Разработанный алгоритм численной оценки распределения веро- 
ятностей для сложных пуассоновских потоков показывает, что в случае 
большого разброса потребностей покупателей в товаре размеры ежеднев- 
ных продаж могут иметь существенно негауссовское распределение. 
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